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ABSTRACT 

Observations of deuterium fractionation in the solar system, and in interstellar and circumstellar 
material, are commonly used to constrain the formation environment of volatiles. Toward protoplan¬ 
etary disks, this approach has been limited by the small number of detected deuterated molecules, 
i.e. DCO+ and DCN. Based on ALMA Cycle 2 observations toward the disk around the T Tauri 
star AS 209, we report the first detection of N 2 D+ (J=3-2) in a protoplanetary disk. These data 
are used together with previous Submillimeter Array observations of N 2 H+ (J=3-2) to estimate a 
disk-averaged D/H ratio of 0.3-0.5, an order of magnitude higher than disk-averaged ratios previously 
derived for DCN/HCN and DCO+/HCO+ around other young stars. The high fractionation in N 2 H+ 
is consistent with model predictions. The presence of abundant N 2 D+ toward AS 209 also suggests 
that N 2 D+ and the N 2 D+/N 2 H+ ratio can be developed into effective probes of deuterium chemistry, 
kinematics, and ionization processes outside the CO snowline of disks. 

Subject headings: astrochemistry—ISM: molecules—protoplanetary disks 


I. INTRODUCTION 

The gas and dust disks surrounding young stars are the 
formation sites of planets. Obtaining a chemical inven¬ 
tory of these disks is key to developing effective molecular 
probes of disk ionization, temperature, and other charac¬ 
teristics that set planet formation efficiencies. Applying 
such probes to disk studies can illuminate the processes 
that led to the present-day composition of solar system 
bodies and exoplanets. 

The thermal histories of solar system, proto¬ 
planetary disk and interstellar material are often 
probed via measurements of deuterated isotopo l ognes 
of common molecules (e. g. lAikawa fc Herbs^ 1 1 9991: 
iMnmma fc CharnlevI 1201 ll) . Molecular D/H ratios in 
cold interstellar and circumstellar environments can be 
several orders of magnitude greater than the co s mic D/H 
abundance of ^ 10 “^ (e.g., iCrapsi et al.ll2QQ5l: iQi et al.l 
[ 200 a lOberg et al.l [20ll iTeagne et al] [201^ . The in¬ 
crease in the D/H ratio of specific molecules, i.e. the 
deuterium fractionation, is a consequence of the lower 
zero-point energies of deuterated molecules compared to 
their non-deuterated isotopolognes. The resulting small 
energy barrier for a molecule to exchange a deuterium 
atom for a hydrogen atom favors deuterium fractionation 
in cold environments. High D/H ratios in the ISM help 
to identify regions that are currently cold, e.g. pre-stellar 
cores, and to infer a low-temperature history for material 
in warmer regions. In the solar system, measurements of 
D/H ratios in comets set constraints upon the tempera¬ 
ture (and therefore th e location) at which they formed in 
the S olar nebula (e.g. iHartogh et al.|[ 2 QTTI:IAltwegg et al.l 
iflii), and the high D/H ratio of terrestrial water indi- 
cates inheritance from the cold parent molecular cloud 
(jCleeves et al.ll^Idbl) . 

Although more than 30 deuterated molecules have 
been detected throughout va rious pre-stella r and 
protostellar environm ents (see iRpueff fc GerinI I2QQ3L 
iCeccarelli et al.ir2QQ7l and references therein), only two of 
them, DCO+ and DCN, had been detected in protoplan¬ 


etary disks (jvan Dishoeck et al.ll2QQ3l: IQi et al.ll 2 QQ 8 f ). A 
third deuterated molecule t hat has been p redicted to be 
abundant in disk s is N 9 D+ (|Willacyll2QQ7l) . N 2 D+ forms 
via the reaction (jPalgarno fc Lepplll984l ) 

H 2 D++ N 2 ^H 2 + N 2 D+. (I) 

The formation of H 2 D+ co mpared to Hj~ is enhanced at 
temp eratures below 50 K (jHerbstI [19821: ICeccarelli et al.l 
120141 ). This leads to an increase in the N 2 D+/N 2 H+ ra¬ 
tio in cold environments. Because H 2 D+ is efficiently 
destroyed by gas-phase CO, abundant N 2 D+, co-existing 
with its non-deuterated isotopologue N 2 H+, is only ex¬ 
pected ou tside the CO sn owline, where CO freezes out 
on grains (jOi et al.ll2013b[) . Together with N 2 H+, N 2 D+ 
may be the best and perhaps only tracer of the chem¬ 
istry and ionization of the midplane in the outer disk re¬ 
gion (jCleeves et al.ll^ldal ). the proposed formation site 
of Kniper Belt objects, comets, Ur anus, and Neptune 
(lOi et al.ll 201 ^ : lAli-Dib et all[ 2 ^ . 

As part of an ALMA Cycle 2 disk survey of deuterated 
molecules, we searched for N 2 D+ toward AS 209, a 1.6 
Myr old K5 T Tauri star in p Ophiuchi located 13 Qd: 50 
pc away (iHerbig fc Belli IM I van Leeuwei] 120071 ). AS 
209 is surrounded by a large and gas-rich disk, extend- 
ing out to 340 AU in CO emission (j Andrews et al.l[ 200 ^ 
iQberg et al.ll 201 l[) . Molecules previously detected in the 
disk include CO, H CO+, DCO+, N 2 H+, HCN, and CN 
(jOberg et al.lDOlTI ). In this letter, we present the first 
detection of N 2 D+ in a protoplanetary disk. We use ex¬ 
isting Submillimeter Array observations of N 2 H+ in the 
same disk to calculate its deuterium fractionation and 
discuss the implications for future disk studies. 

2. OBSERVATIONS 

AS 209 (J2000.0 R.A. 16N9“15®.29, decl. 

-14°22'08".6) was observed 2014 July 2 during 
ALMA Cycle 2 with 21 minutes of on-source integration 
time. Thirty-four 12 m antennae were used for the 
observations, with a maximum baseline of 650.3 m. 
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Thirteen spectral windows (SPWs) were observed in 
Band 6, each with 61.04 kHz resolution. The N 2 D+ 
J=3-2 line was in SPW 9, centered on 231.32183 GHz, 
and the J=2-l line (used in this study as a ylsr 

and spectral profile reference) was in SPW 7, centered 
on 220.39868 GHz. The 1.4 mm dust continuum flux 
was extracted by averaging the channels in five line-free 
SPWs. Each of these windows was 59 MHz wide. 

ALMA/NAASG staff performed bandpass and phase 
calibration with the quasar J1733-1304, and the flux was 
calibrated with Titan. The uncertainty in the absolute 
flux is ^10%. The calibrated visibilities were decon¬ 
volved and GLEANed with the GASA software package 
(version 4.2.2). The integrated flux density of the 1.4 
mm dust continuum is 242T24 mJy. This was measured 
by integrating over the area of the disk delimited by a 
2cr contour, where cr=0.3 mJy beam“^ was the contin¬ 
uum rms measured in a signal-free region of the intensity 
map. After continuum subtraction in the uv-plane, the 
data cubes for the SPWs were rebinned to 0.4 km/s. 
To increase signal to noise, the visibilities for the N 2 D+ 
J=3-2 line were tapered to 1". Because the spatial distri¬ 
bution of the emission changes between different velocity 
channels, GLEANing for both N 2 D+ and was per¬ 

formed using a Keplerian rotation mask based on the 
emission patterns of the stronger line, ^^GO 2-1. The 
channel rms for ^^GO 2-1 was 17 mJy beam“^, while the 
rms for N 2 D+ 3-2 was 8.5 mJy beam“^. 

To estimate a D/H ratio for N 2 D+/ N 2 H+, this pa¬ 
per also uses AS 209 data originally published in Oberg 
et al.’ s Submillimeter Array protoplanetary disk survey 
The 1.4 mm dust continuum flux density from 
their dataset (also obtained by integrating over the area 
of the disk encompassed by the 2cr contour) is 219 ± 
22 mJy. Thus, within the uncertainties, the flux density 
from the ALMA observations is consistent with that of 
the SMA observations. 
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Figure 1. Continuum and line emission maps toward AS 209. 
The topmost panel shows the 1.4 mm dust continuum emission, 
with contours at [5, 10, 20, 40, 80, 160] <7, where a = 0.3 mJy 
beam“^. Integrated intensity maps (summed from 2.6 to 7.0 km 
s“^) with contours at [2, 3, 4, 8] <7 are shown in the left column 
for 2-1 , with a = 29 mJy km s“^ beam“^, and N 2 D+ 3-2, 

with a = 12 mJy km s“^ beam“^. The corresponding panels in the 
right column show the integrated intensities for the blueshifted (2.6 
to 4.6 km s“^) and redshifted (5.0 to 7.0 km s“^) emission. The 
pointing center in each panel is marked by a cross, with position 
offsets in arcseconds. Synthesized beams are drawn in the lower 
left of each panel. 


3. RESULTS 

Eigure[T] shows the AS 209 1.4 mm dust continuum, in¬ 
tegrated intensity maps (spanning channels from 2.6 to 
7.0 km s“^), and blue and redshifted disk emission maps 
for the N 2 D^ 3-2 and 2-1 lines. The CO emission 

is symmetric, whereas most of the emission above the 2cr 
level for N 2 D+ 3-2 is east of the phase center. The peak 
intensity of N 2 D+ is 3.3cr, where a=12 mJy beam“^ km 
s“^ is the rms measured in a signal-free portion of the 
integrated intensity map. The velocity field of ^^CO 2-1 
in Eig. [T] is well-resolved and consistent with the Kep¬ 
lerian rotation of a disk. The redshifted N 2 D+ emission 
is consistent with the redshifted emission, but no 

blueshifted emission is observed for N 2 D+ above the 2cr 
level. This is spatially consistent with the asymmetry 
observed in the total integrated intensity map. While 
the morphologies of the contours differ between the to¬ 
tal integrated intensity maps and the redshifted emission 
maps for N 2 D+ due to the differing noise contributions 
for the different velocity ranges, the mean intensity mea¬ 
sured from the integrated intensity map is within 5% of 
the mean intensity measured over the redshifted channels 
only, confirming that most of the emission is redshifted. 

Velocity channel maps for N 2 D+ 3-2 are shown in Eig. 
[21 Line fiux at >3cr level is detected near the phase center 
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Figure 2. Velocity channel maps of N 2 D+ 3-2 emission toward 
AS 209, binned by 0.4 km s“^. 2<7 and 3<7 contours are drawn, 
with <7 = 8.5 mJy beam“^. Channel velocities [km s“^] are listed 
in the lower right corner of each panel. The pointing center is 
denoted by a cross, with offsets from the pointing center indicated 
in arcseconds on the axes. Synthesized beams are drawn in the 
lower left of each panel. 
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in three channels, and at the 2cr level in another three 
channels. The progression of emission from west to east 
of the pointing center is consistent with disk rotation, 
confirming the disk origin of the N 2 D+ emission shown 
in Fig. [H We note that although blueshifted emission is 
visible in the velocity channel maps, it is not apparent in 
the integrated intensity map in Fig. [T] above the 2cr level 
because the map is integrated across a velocity range 
that includes a number of channels without blueshifted 
emission. 

Figure [ 3 ] displays spectra for the N 2 D+ 3-2 and 
2-1 lines, extracted using the Keplerian CLEAN masks, 
together with the N 2 H+ 3-2 SMA spectrum. The double- 
peaked profile of 2-1 is typical of a rotating disk. 

The spectrum of N 2 H+ 3-2, though with a lower signal- 
to-noise ratio (SNR), is also consistent with a double- 
peaked profile. As with the integrated intensity maps, 
the emission for N 2 D+ is asymmetric, with only the red- 
shifted peak clearly visible. However, the position of the 
N 2 D+ peak is consistent with the positions of the red- 
shifted peaks for N 2 H+ and ^^CO. The hyperfine struc¬ 
ture is not resolved for N 2 D+. 

The integrated fluxes are given in Table [H The un¬ 
certainty on the integrated flux was calculated using the 
formula a = rms x FWHM/ where the rms is mea¬ 

sured from a featureless portion of the spectrum for each 
transition, the FWHM is measured from the 2-1 

line, and rich is the number of channels spanning the 
FWHM. 

Using the N 2 H+ line data from lOberg et al.l (|2Qllf ) 
in conjunction with the ALMA observations of N 2 D+, 
we estimate a disk-averaged N 2 D+/N 2 H+ ratio toward 
AS 209. Because the SNR is low for both molecules, 
their spatial distributions are not well-constrained by the 
current set of observations. Hence we assume that the 
molecules are co-spatial in the disk based on their sim¬ 
ilar formation and destruction chemistry. In particular, 
the N 2 D+ and N 2 H+ distributions are both expected to 
be regulated by destructive reactions with CO gas, and 
therefore the degree of freeze-out of CO and N 2 . The 
N 2 H+ intensity, averaged over a 5" x 5" circle centered 
on the continuum peak, is 320T40 mJy beam“^ km s“^, 
while the N 2 D+ intensity, averaged over a 4" x 4" circle 
centered on the continuum peak, is 9±2 mJy beam“^ km 
s“^. In each case, the surface area over which the inten¬ 
sity was averaged was based on the extent of the molec¬ 
ular emission above the 2 x rms level and was largely set 
by the synthesized beam size. 

Column densities were derived from the integrated 
i ntens ities using the LTE formula from iRemiian et al.l 
(I 2001 . in the case where the emission is optically thin 
and the continuum brightness temperature is assumed to 
be much lower than the excitation temperature: 


N^ 2.04fAI,.Q 


OaOb 




( 2 ) 


j A/ dv ( Jy beam“^ km s“^) is the integrated line 
intensity, Oa and 0^ (arcsec) are the major and minor 
axes of the ellipse defined by the EWHM of the Gaus¬ 
sian beam, and Tq (K) is the excitation temperature. 
The molecular parameters, which are the rotational par¬ 
tition function Qrot, transition upper state energy E^^ 
(K), transition rest frequency u (GHz), line strength S', 


and electric dipole moment /i (Debye), were all obtained 
from the Cologne Datab ase for Molecular Spectroscopy 

(|Mii]]CT et al.ll2^l2nn^ . 

The excitation temperature is expected to be between 
10 and 25 K in the outer regions of the d isk where N 2 D^ 
and N 2 H+ reside (e.g. IWillacvl I2QQ7I: iCeccarelli et al.l 
l2Q14f ). At 10 K, the disk-averaged column density of 
N 2 D+ is 2.1 ± 0.5 X 10^^ cm“^ and that of N 2 H+ is 
6.3 ± 0.8 X 10^^ cm“^. At 25 K, the disk-averaged col¬ 
umn density of N 2 D+ is 1.4±0.3 x 10^^ cm“^ and that of 
N 2 H+ is 3.1 ± 0.4 X 10^^ cm“^. At a given temperature, 
the dominant source of uncertainty in the LTE approxi¬ 
mation for the column density is due to the uncertainty 
on the intensity measurement. The resulting D/H ratio 
is estimated to be between 0.3 and 0.5. As long as the 
N 2 D+ and N 2 H+ are co-spatial, the ratios of their disk- 
averaged column densities can be estimated with greater 
confidence than their absolute column densities, which 
are more sensitive to the assumed surface area of the 
emission. 
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Figure 3. Spectra of the 2-1 and N 2 D+ 3-2 transitions (ex¬ 

tracted with Keplerian masks) toward AS 209, as well as a spec¬ 
trum of N 2 H+ 3-2 toward AS 209 from Submillimeter Array data 
published in lOberg et a DCoU). The vlsr is marked by the ver¬ 
tical dashed lines. 


To check the validity of the optically thin LTE approx¬ 
imation for estimating the column densities of N 2 H+ and 
N 2 D+ toward AS 209, column densities for N 2 H+ were 
also calculated with RADEX, a ID no n-LTE radiative 
transfer code (|van der Tak et af] I2QQ7I ). To cover the 
range of typical disk conditions, we ran a grid of mod¬ 
els with kinetic temperatures between 10 and 25 K and 
an H 2 number density (representing the dominant colli¬ 
sion partner) of 10^-10^^ cm“^. The resulting column 
densities ranged from 1.7 to 7.8x10^^ cm“^, similar to 
the column density range estimated through LTE calcu¬ 
lations. 


4. DISCUSSION 

4.1. N 2 H^ deuterium fraetionation 
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Table 1 

Summary of Line Observations 


Transition 

Rest Frequency 
(GHz) 

(K) 

Beam 

Integrated Flux 
( mJy km s“^) 

Channel rms 
( mJy beam“^ ) 

i3CO 3=2-1^ 

220.39868 

15.9 

l."02x 0."74 (87°.18) 

1940[37] 

17 

N 2 D+ 3=3-2^ 

231.32183 

22.2 

l."05x 0."99 <°.83) 

80[18] 

8.5 

N 2 H+ 3 = 3 - 2 '^ 

279.51173 

26.8 

3."74x 2."50 (19°.26) 

597(78] 

80 


^ From ALMA observations. 

^ From data originally published in lOberg et al.l (|2011IL 


The estimated N 2 D+/N 2 H+ disk-averaged abundance 
ratio of 0.3-0. 5 is consistent with the disk models of 
iWillacvl (1200^ . in which the D/H ratio was predicted 
to range between 0.35 and 7.0 in the outer disk. The 
predicted range reflects the dependence of deuterium 
chemistry on desorption from dust grains due to cos¬ 
mic rays, thermal proce sses, and UV r adiation, which 
varied between models. IWillacvl (j2007[ ) also predicted 
that the D/H ratio for N 2 D+ in the outer disk would be 
one or two orders of magnitude greater than the cor¬ 
responding ratios for DCN and DCO+. This is con¬ 
sistent with our finding that the disk-averaged D/H 
ratio for N 2 H+ toward AS 209 is at least an order 
of magnitude higher than the disk-averaged ratios de¬ 
rived for HCN and HCO+ toward other disks. A disk- 
averaged DCN/HCN abundance ratio of 0.017 was de¬ 
rived for TW Hya (jOberg et al.l l2012l) . Disk-averaged 
DCO+/HCO+ abundance ratios have been estimated to 
be 4x10“^ toward DM Tau, 0.02 toward HD 163296 , 
and 0.035 toward TW Hva (Ivan Dishoeck et al.l 120031 : 
iGuilloteau et al.l 120061: iMathews et al.ll2013f ). However, 
analyses of spatially resolved data for these three disks 
have found local DCO^/HCO^ ratios rising as high as 
0.1-0.3, with radial variations in the a bundances due 
to radial changes in the temperature (|Qi et al.l 120081 : 
IMathews et al.ll20T^ iTeague et al.ll201^ . The local en¬ 
hancement levels revealed by the spatially resolved data 
indicates the value of obtaining spatially resolved data 
for other deuterated lines, such as N 2 D+. 

The high D/H ratio for N 2 D+/N 2 H+ toward AS 209 

and the lower DCO+/HCO^ ratios found in other pro¬ 
toplanetary disks suggest that the deuterium chemistry 
in the outer disk resembles that of dense cores and 
protostellar envelopes. A survey of low-mass starless 
core s found N 2 D^/N 2 H+ ratios ranging from <.02 to 
.44 (|CraDsi et al.l 1200^ . and a survey of Class 0 proto¬ 
stars ranging from <.034 to 

.27 (|EmDrechtinger et al.l 120091) . For the Class 0 proto¬ 
stars in which both N 2 D+ and DCO^ have been mea¬ 
sured, the N 2 D+/N 2 H+ ratio often exc eeds the DCO+/ 
HCO^ ratio by an order of mag nitude (j Jorgensen et al.l 
|2004 lEmprechtinger et all 120091) . S imilarly , toward the 
dense core L1544, Caselli et al. (j2002br ) found that 
N 2 D+/N 2 H+ was ^0.2, versus a ratio of 0.04 for DCO+/ 
HCO+. 

Civen that models posit similar formation pathways 
for N 2 D+ and DCO+, it may be surprising that such 
different D/H ratios are predicted and observed for the 
two molecules. This difference arises because gas-phase 
CO, a parent molecule of DCO+, destroys H 2 D+one_of 
the main parent molecules of N 2 H+ and N 2 D+ (jOi et al.l 


l2013al) . Due to this destruction pathway, N 2 D+ is only 
abundant at low temperatures, where deuterium frac¬ 
tionation is especially efficient and gas-phase CO is de¬ 
pleted due to freezeout onto grains. 

4.2. A 2 D+ emission asymmetry 

The N 2 D+ emission toward AS 209 appears to be 
asymmetric. Higher SNR observations are required 
to confirm the reality of this asymmetry, but it is 
worth considering possible causes for this tentative 
azimuthal N 2 D+ variation. Previously observed az¬ 
imuthal variations in disk molecular emission have 
been ascribed to foreground absorption, projection ef¬ 
fects due to disk inclination, and dust traps cooling 
nearby gas fe.g. lOberg et al.ll2Qlll : lRpsenfeld et al.ll2Q13l: 

Ivan der Plas et M J 1211141). Foreground absorption is un¬ 
likely to cause asymmetric emission for N 2 D+ due to 
the molecule’s low abundance and high critical density. 
For molecular emission originating close to the cold mid¬ 
plane, as is expected for N 2 D^, projection effec ts should 
also play a minimal role (jRpsenfeld et al.ll2Q13f ). Asym¬ 
metry in the dust structure is an unlikely cause, given 
the axisymmetric appearance of the 1.4 mm continuum 
emission shown in Fig. [TJ However, because the N 2 D+ 
distribution is expected to be especially sensitive to ion¬ 
ization and temperature, possible causes for asymmetry 
include X-ray flares or azimuthal density gradients in 
the o uter disk (e.g. iBirnstiel et al.l I2Q13I: iCleeves et al.l 
I2Q15D . Temperature asymmetries could be probed with 
D/H ratios derived from high SNR spatially-resolved 
observations of N 2 D+ and N 2 H+, while asymmetry in 
the midplane ionization structure could be evaluated by 
comparing N 2 D+ to neutral tracers such as H 2 CO, an- 
other molecule e xpected to form past the CO snowline 
(IQi et al.ll20T^ . 

4.3. Developing N 2 D^ as a disk probe 

The similarity between the deuterium fractionation in 
disks and interstellar cloud cores suggest that analogous 
to its use as a tracer of CO-depleted cold gas in cores 
(|Caselli et al.ll2QQ^ lEmprechtinger et al.|[2QQ9D . N 2 D+ 
with N 2 H+ could be developed into a unique tracer of 
the cold, CO-depleted outer disk midplane. In addition, 
measurements of N 2 D+ and N 2 H+, together with DCO+ 
and HCO+, have been used to estimate the ion ization 
fraction of cloud cores (|Caselli et al.l 1199^. l2QQ2bl ). This 
approach could be adapted to evaluate the ionization 
fraction of disks. 

Spatially resolved observations through ALMA would 
allow the D/H ratio as a function of radius to be derived, 
thereby providing a stronger test for chemical theories re¬ 
garding how N 2 D+ is formed and destroyed throughout 
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disks. A sensitive survey of N 2 D+ and N 2 H+ in disks 
would also be useful to determine the relationship be¬ 
tween deuterium fractionation and other properties of 
the disk, such as mass, age, and ionization processes. 

As molecules with abundances anti-correlated with the 
abundances of CO isotopologues and related species (e.g. 
HCO+, DCO+), N 2 D+ and N 2 H+ complement the cur¬ 
rent set of common disk tracers. While previously de¬ 
tected disk molecules have been abundant in the in¬ 
ner disk and in the warm molecular layer high above 
the midplane, N 2 D+ and N 2 H+ are suitable for tracing 
kinematics, ionization, and temperature structure in the 
outer disk midplane. Although H 2 D+ is also expected to 
be abundant in the outer disk midplane, it has not yet 
been detected in disks and is expecte d to require very 
long integration (jCleeves et al.l [2QT4al ). Therefore, fur¬ 
ther observations and modeling of N 2 D+ and N 2 H+ are 
key to developing a fuller picture of how disk character¬ 
istics vary spatially. Constraining such properties can 
clarify how volatiles have been formed, transported, or 
destroyed in disks, giving rise to the planets we observe 
now. 

This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2013.1.00226. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRC (Canada) 
and NSC and ASIAA (Taiwan), in cooperation with the 
Republic of Chile. The Joint ALMA Observatory is oper¬ 
ated by ESO, AUI/NRAO and NAOJ. The National Ra¬ 
dio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. We thank Adam 
Leroy and the NAASC for assistance with calibration and 
imaging, Ryan Loomis and Viviana Guzman for helpful 
discussions on the data, and the referee for useful com¬ 
ments on the paper. K. 1. O. also acknowledges funding 
from the Simons Collaboration on the Origins of Life 
(SCOL), the Alfred P. Sloan Foundation, and the David 
and Lucile Packard Foundation. 
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